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Abstract 

The generalization of the Zubarev nonequihbrium statistical operator (NSO) method for the case 
of Renyi statistics is proposed when the relevant statistical operator (or distribution function) is 
obtained based on the principle of maximum for the Renyi entropy. The nonequilibrium statistical 
operator and corresponding generalized transport equations for the reduced-description parameters 
are obtained. A consistent description of kinetic and hydrodynamic processes in the system of 
interacting particles is considered as an example. 
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I. INTRODUCTION 



In investigations of complex self-organizing, fractal structures and various physical phe- 
nomena such as subdiffusion, turbulence, chemical reactions as well as various economical, 
social and biological systems Gibbs distribution function does not provide agreement with 
observable phenomena. For these systems the power distributions are inherent jl}]. They 
can not be obtained from the maximum entropy principle for the Boltzmann- Gibbs entropy 
underlaying both an equilibrium and nonequilibrium statistical thermodynamics {2 4|. 

In papers of A.G. Bashkirov js-Sj the use of the Renyi [o], [loi entropy as statistical 
entropy for investigation of complex systems is proposed. It depends on the parameter q 
(0 < g < 1) and at g = 1 it coincides with Boltzmann-Gibbs entropy. Based on the maximum 
entropy principle for the Renyi entropy the power Renyi distribution was obtained in the 
case of equilibrium. At g = 1 it develops into the Gibbs canonical distribution. Herewith 
?7 = 1 — g is considered as an order parameter. With the increase of it the statistical 
Renyi entropy grows to its maximum which the power Renyi distribution corresponds to. 
Moreover the Renyi entropy derivative with respect to 7] suffers a sudden change at r/ = 0, 
that is a sort of phase transition to more ordered equilibrium state takes place. The papers 



of Abe 
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12|, Arimitsu [l3|,l]J, Lesche jl5[, Masi [l6| and others were devoted to axiomatic 



argumentation and problems of stability of the Renyi entropy and its linearized form — the 
Havrda-Charvat-Tsallis entropy |l7. 



in books 



18l |. The properties of the Renyi entropy are discussed 



as well. Nowadays the Tsallis entropy is widely applied in various 



directions of nonextensive statistical mechanics 
of subdiffusion [2^ and turbulence 
gases and plasma [27 



25 



2l|-23|. The examples are the phenomena 



261], the investigations of transport coefficients in 
28l ] as well as quantum dissipative systems The problems of 



the construction of equilibrium thermodynamics in the framework of the Renyi and Tsallis 



statistics were discussed in 



30 



-|37|. The energy fluctuations 33|], kinetics of nonequilibrium 



ozone layer [40|, universality in non- 
431 were investigated within the Tsallis 



plasma [381], problems of self-gravitating systems 
Debye relaxation 4l| and complex systems 42, 
statistics. 

In the present paper one approach of formulation of an extensive statistical mechanics 
of nonequilibrium processes is considered based on the nonequilibrium statistical operator 
method by Zubarev and the maximum entropy principle for the Renyi entropy. The 
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consistent description of kinetic and hydrodynamic processes in the system of classical in- 
teracting particles is considered as an example. 



II. MAXIMUM ENTROPY PRINCIPLE FOR THE RENYI ENTROPY AND 
NONEQUILIBRIUM STATISTICAL OPERATOR 

Nonequilibrium state of a classical or quantum system of interacting particles is com- 
pletely described by the nonequilibrium statistical operator (nonequilibrium distribution 
function) g{x^;t). The latter satisfies the classical or quantum Liouville equation 



^g{x'';t)+tLNg{x'';t) = 0. (1) 

iLjs! is the Liouville operator the system of interacting particles which in classical case has 
the following form: 

a,^±EL.^-i±^^r„)(±-^). (2) 

dri 2 i^^dri \dpi dpj J 

Here xj = {pj,rj} is the phase variables of j-particle, ^{rij) is the interaction energy of two 
particles, pj the j-particle momentum and m its mass, rij = \ri — fj\ the distance between 
pair of interacting particles. 

The function q[x^ ; t) is symmetric regarding to inversion of phase variables of any pair 
of particles xi ^ xj and satisfies the normalization condition J drNg{x'^;t) = 1, dT = 
(dx)^ /N\, dx = dpdr. 

In order to solve the Liouville equation ([T]) we will use the Zubarev nonequilibrium sta- 
tistical operator method 2rM]- Within its framework we will be looking for solutions of the 
equation ([T]), which are independent of the initial conditions. The solutions will depend on 
time explicitly only, i.e. through the observable quantities 

dTr,PnQ{x''-t) = {P^Y (3) 

selected for the reduced description of nonequilibrium states of the system. In particular, for 
the description of the hydrodynamic state of the system of classical interacting particles the 
nonequilibrium averaged values of densities of number of particles (n(r*))*, their momentum 
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(i){f)y and total energy (e(r))* can be chosen as the parameters of the reduced description 
^. They satisfy the respective conservation laws. In the case of kinetic description the one- 
particle fi{p,r;t) and two-particle f2{p,r,'p',r';t) distribution functions can serve as the 
reduced-description parameters. For the investigation of properties of magnetic and polar 
systems the averaged values of densities of magnetic (m(r))* and dipole {d{T)Y moments 
can be used respectively besides the hydrodynamic variables. 

When the basic parameters of the reduced description are defined, the solution ^(x^; t) 
can be presented by means of the nonequilibrium statistical operator method in the following 



general form 



2i-y 



Qix^'-t) = Qrelix^'-.t) - / c'^'' ''^T {t,t') (1 - Prel{t')) iL^ Qreli^'' ■,t')dt' , (4) 



oo 



where T{t,t') = exp_|_|— J^* (1 — Preiit')) iL^dt'^ is the evolution operator containing pro- 
jection; exp_,_ the ordered exponential. Preiit') is the generalized Kawasaki-Gunton projection 
operator whose structure depends on the form of the relevant statistical operator (distribu- 
tion function) grei{x^]t). The latter will be determined using maximum entropy principle 
for the Renyi entropy 

S''{g) = -^\n [ dT^g^t). (5) 



1-q 

The corresponding functional at fixed parameters of the reduced description with taking 
into account the normalization condition has the form 



L'^ig) = Y^^ln j dT^g^it) -a j dTNg{t) -J^F^Xt) j dV^Pngit) 



(6) 



Fn{t) are the Lagrange multipliers. Equalizing its functional derivative to zero we obtain 
the relevant statistical operator corresponding to the Renyi entropy maximum: 



Qrelit) = ^ - ^ E Fn{t)6Pnj , (7) 

Znit) = j dV^ ^1 - ^ ^ F„(t)5P„ j ' . (8) 

Zji(t) is the partition function of the relevant statistical operator, 6Pn = Pn — {PnY, and 
the parameter a in ([6]) is determined by the relation 

a = -^^-J2^nmPnY- (9) 



q-l 



The Lagrange multipliers in ©-(E]) are defined from the self-consistency conditions: 



(10) 



(■ ■ -Yrei — J dr^ . . . Qrei{x^ ] t) . Since the relevant statistical operator is now known for the 
basic set of the reduced-description parameters, we can obtain the nonequilibrium statistical 
operator, establishing the structure of the projection operator: 



Prel{t)g' = Qrelit) - ^ . 



(11) 



The variation derivative of the relevant statistical operator in (|TT1) can be presented in the 



form: 



S Qrelit) 
S{PmY 



Qrel{t)5 



HPmY 



5P„. 



(12) 



where 



.] = [.. .]-([...]) 



reZ' 



and we use the notation 



m = 1 
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Y,Fn{t)5Pn. 



(13) 



(14) 



The Lagrange multipliers derivatives with regard to the reduced-description parameters 
6Fn(t)/S{PmY we calculate in the following way 

-1 

(15) 



5F„(t) / S{PmY 



It can be done in general case. Thus 

SJPrnY_SJP^_ f ^ 

After calculating 6grei{t)/SFn{t) in the right-hand side of the relation flTB]) . we obtain the 
set of equations for desired derivatives 



5{PmY _ 1 

5F„(t) q 



^ SF^{t) q 



I reV 



(17) 



Its solution in the matrix form is 



5{py 



I - -{6Pr\t))i,F 



5F(t) 

where / is the unit matrix and 

5F„(i) [SF{t) 
Thus the functional derivative can be written in the form: 

Sgrel{t) 



-{5Pr\t)5P)\,i = fit), 



= Qrel{t)5 



Then the Kawasaki-Gunton projection operator has the following structure: 

Prel{t)g' = Qrel{t) J dTNlg'} 



(18) 



(19) 



(20) 



(21) 



+ ^Qrel{t)S 



-r\t)lFUt)+J2f^nit)SPn 



X 



J dTNiPruQ'} - {PrnY j dTn{Q'] 



It is further necessary to explore an action of the operators Prei{t)iL]si on the relevant 
statistical operator. Since 
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then 



iLNQrel{t) = -Qrel{t)-^ ^{t)^Fn{t)Pn = A{t) Qrel{t) ■, 



Prel{t)iLNQrel{t) = P {t) A{t) Qrel{t) = / dV {A{t) Qrel{t)} 



(22) 



(23) 



-r\t) (F^(i) + $^/-^(i)5P„ 



X / dTM{PmA{t)Qrel{t)} - {PmY / dV N{A{t) Qrel{t)} , 



where 



j drN{PmA{t)grel{t)} - {PmY j dV N{A{t) Qrelit)} = {5PmA{t)Yrel- (24) 
Thus Prel{t)iLNQrel{t) = Prel{t) A{t) Qrel{t) = {P (t) A{t)) Qrelit) , whcrC P{t) is the projcction 

operator which now acts on dynamic variables: 



Pit) ■■■=(■■ -Yrel 



(25) 



V^(t) (Fm{t)+Y,f^Ut)SPn 



(• • • SPjnYrel- 



So far as 



A(t) = — V>-i(t)VF„(t)P„ 
q ^-^ 



we have 



p{t)A{t) = — Y,Fr,{t)i,r\t)Pn) 



t 

rel 



(26) 



(27) 



X 



- VF„(t)(^-l(t)P„) J (P,n - {Pm) 

q ^-^ J 



rel ) 



rel 



Considering (l22|) -(l27j) we can present (1 — Preiit))iLN Qreiit) in the form: 

(1 - Prelit))iLN0rel{t) = 

= (1 - P{t))lLNgrel{t) = -^In{t)Fn{t)Qrel{t), 



(28) 



where 



Ut) = {l-P{t))-i;-\t)K 



(29) 



are the generahzed flows. Taking into account (128|) we can now write down an exphcit 
expression for the nonequihbrium statistical operator 



g{x ;t) = Qreiix ;t) 



(30) 



f e^^*'-*)T(^,^')/n(^')^n(^')^? 



rel{x''-t')dt' 



It allows us to obtain the generalized transport equations for the reduced-description 
parameters. They can be presented in the form: 

^(P^)*= + e^^*'~*Vmn(t,t')^n(Orf^', (31) 

with the generalized transport kernels (memory functions) 



{t,t')= / dTN{PmT{t,t')Irr{t')Qrel{t')} 



(32) 



which describe the dissipative processes in the system. 



III. GENERALIZED TRANSPORT EQUATIONS FOR A CONSISTENT DE- 
SCRIPTION OF KINETICS AND HYDRODYNAMICS IN THE RENYI STATIS- 
TICS 



For a consistent description of kinetic and hydrodynamic processes in classical (or quan- 
tum) systems of particles interacting in the volume V the nonequilibrium one-particle 
distribution function fi{x]t) = (ni(a;))* and averaged value of potential energy of inter- 
action eint{r;t) = (emt(r))* are the basic parameters of the reduced description The 
last is defined through the two-particle distribution function f2{x,x';t) = (n2(x, x'))*. Here 
ni(x) and h2{x,x') are the phase densities of the microscopic distribution of particles, and 
^intif) = ^fdpfdp'f dr'^{\r — r'\)h2{x,x'). In this case according to ([7]) we obtain the 
relevant statistical operator 



1 / q-1 
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Qrei{t) = — yl q~\J + J dxa{x;t)5ni{x;t)jj , (33) 

where 

ZR{t) = J dTN (^l - ^-^{ j dfp{f;t)6iint{f;t) + j dxa{x-t)5hi{x]t)^^ " ' (34) 

is the partition function. The parameters P{r;t) and a{x;t) are defined from the self- 
consistency conditions: 

According to (jl]) and taking into consideration ( 133|) we write down the nonequilibrium 
statistical operator of a consistent description of kinetic and hydrodynamic processes in the 
system: 

gix""; t) = greiix""; t) + f e'^''~'^T{t, t') (36) 

J — oo 

X (^j dr^l3{f-t')r;'\r-t') + j dx'a{x';t')Ir,{x';t')^ Qrei{x^ ]t')dt' . 

Here 

In{x'-t') = (1 - P{t'))Kp-\t')iLNhi{x') (37) 
IT^f'-t') = (1 - P{t'))-r\t')iLj,e,r^t{n (38) 
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are the generalized flows. With the help of the nonequilibrium statistical operator obtained 
one can derive the generalized transport equations for the basic set of the reduced- description 
parameters according to (pii) : 

^ (ni(x))* = - [ dx'(^nnix,x';t)a{x';t) (39) 



+ ^ y df^ne{x, r; t)/3{f; + j dx' j e^(*'"*Vnn(a;, t, t')a{x'] t')dt' 
+ ! df' f e'^''-'^^U^,f';t,t')l3{r';t')dt', 



oo 



d 1 

{£intir)y = - I dx'^eh{r,x';t)a{x';t) (40) 



Ot q 

+ IJ dr'^eeir, r'- t)P{f- 1) + J dx' J e"(*'-*Ven(r, x'; t, t')a{x'; t')dt' 
+ ! rfr" f e<''-'^<fUry;t,t')f3{r';t')dt'. 



The functions 



^nh{x,x';t) = J dTNni{x)ilj'^{t)iLNni{x')grei{x'^;t), (41) 
^nsix,r;t) = j dTNhi{x)i)~^{t)iLNiint{r)Qrei{x^ ;t), (42) 
$e^(f,x';t) = drNeint{r)ip~^{t)iLNni{x')grei{x^]t), (43) 



(!>ee{r,f^;t) = j dTNeint{r)i) ^{t)iLNeint{^)Qrei{x^ ]t), (44) 

are the time correlation one obtained by means of the relevant distribution and which contain 
the function 

il){t) = 1 — ( / dfp{f;t)6iint{r;t) + / dxa{x;t)6ni{x;t)] , (45) 



At q = 1 tplt) = 1, we have the transition to the relevant Gibbs distribution when the 
nonequilibrium one-particle distribution function fi{x;t) = (ni(x))* and averaged value of 
potential energy of interaction eint{r]t) = {eint{r)Y are the parameters of the reduced de- 
scription [44]. The generalized transport kernels ^Pnn{x, x'\ t, t'), fne{x, r'; t, t'), ^enif^ x'\ t, t'), 
iPee{r,r';t,t') have the structure of fl5^ with the corresponding flows ([ST]), f l55]) . 

For the case when the nonequilibrium one-particle distribution function fi{x; t) = (fii(a;))* 
is the single parameter of the reduced description (the contribution of nonequilibrium aver- 
aged potential energy of interaction is much smaller then the kinetic energy, e.g. the case 



of rare gases) the set of transport equations (l39!) - (H0l) reduces to the generahzed kinetic 
equation: 

1 



(x, x'; t)a{x'; t) 



(46) 



+ / dx' 



3^^*'"*Vnn(a;, x'; t, t')a{x'; t')dt' . 



At g = 1 it transforms into the kinetic equation of 4J] with the transport kernel calculated 
using relevant distribution function Qrei{t) = Y^=i In this case at g = 1 within the 

NSO method j2-4| the Liouville equation ([1]) should be solved with the boundary condition 



d_ 

di 



N 



g{x ; t) + iLi^Q^x ; t) 



-e g[x 



N 



(47) 



It corresponds to the Bogolyubov hypothesis of weakening of the correlations be- 
tween particles. Integration of this equation with respect to phase variables J dT^-i, 
J dV is[-2----j dV isi-s leads to the BBGKY hierarchy for the nonequilibrium particle distribu- 
tion functions. The question about the BBGKY hierarchy when in the boundary condition 
( 1i7|l Qreiit) is equal to ( 133|) is interesting. With q = I Qreiif) transforms into the Gibbs form 



4^ which 



and we obtain the BBGKY hierarchy with the modified boundary condition |3, 
takes into account many-particle correlations. It allows one to obtain in the pair collision 
approximation the revised Enskog theory and the Enskog-Landau kinetic equations for the 
neutral and charged hard sphere system, respectively 3I.I4. |44-|46|. 



IV. CONCLUSIONS 



For the nonequilibrium system of interacting particles within the framework of the 
Zubarev NSO method we obtained the nonequilibrium statistical operator Q{t). It satis- 
fies the Liouville equation with the boundary condition describing the relaxation of the 
NSO to the relevant statistical operator Qreiit). The latter is constructed based on the 
maximum entropy principle for the Renyi entropy at fixed values of the reduced-description 
parameters (-Pn)* with taking into account the normalization condition. By means of the 
NSO the generalized transport equations for the parameters of the reduced description {PnY 
are obtained with regard to Kawasaki- Gunton and Mori projection. Such an approach is 
applied to a consistent description of kinetic and hydrodjTiamic processes in the system 
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of classical interacting particles. As a result both the nonequilibrium statistical operator 
and the generalized transport equations are obtained, when the nonequilibrium one-particle 
distribution function fi{x;t) = (ni(a;))* along with the nonequilibrium averaged value of 
the potential energy of interaction eint{r;t) are selected as the reduced-description parame- 
ters. At q = 1 the known results based on the Gibbs statistics are reproduced. Naturally, 
an interesting question about the investigation of time correlation functions and transport 
coefficients based on the NSO method within the Renyi statistics arises. 
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